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Abstract. An image restoration approach based on a Bayesian maximum entropy 
method (MEM) has been applied to a radiological image deconvolution problem, that 
of reduction of geometric blurring in magnification mammography. The aim of the 
work is to demonstrate an improvement in image spatial resolution in realistic noisy 
radiological images with no associated penalty in terms of reduction in the signal- 
to-noise ratio perceived by the observer. Images of the TORMAM mammographic 
image quality phantom were recorded using the standard magnification settings 
of 1.8 magnification/fine focus and also at 1.8 magnification/broad focus and 3.0 
magnification/fine focus; the latter two arrangements would normally give rise to 
unacceptable geometric blurring. Measured point-spread functions were used in 
conjunction with the MEM image processing to de-blur these images. The results 
are presented as comparative images of phantom test features and as observer scores 
for the raw and processed images. Visualization of high resolution features and the 
total image scores for the test phantom were improved by the application of the MEM 
processing. It is argued that this successful demonstration of image de-blurring in 
noisy radiological images offers the possibility of weakening the link between focal spot 
size and geometric blurring in radiology, thus opening up new approaches to system 
optimization. 



1. Introduction 



The design of radiological imaging equipment has developed to reflect the best 
compromise between a number of contradictory performance requirements. In particular 
the design of x-ray tubes has evolved to best fit the requirements for the smallest possible 
focal spot size combined with the maximum target heat capacity. In this case the 
performance compromise is between geometrical blurring introduced by a finite focal 
spot size and movement blurring introduced by the long radiographic exposure time 
required to operate a very small focal spot within its heat loading limits. The most 
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common resulting design, rotating target with an acute target angle, represents an 
impressive engineering solution to this problem. 

An alternative potential solution to the problem of geometrical blurring, applicable 
to modern digital receptors, is digital image processing. Provided the characteristics 
of the image blurring function are known, an image can, in theory, be restored to 
that which would have been produced by a perfect point focal spot. However, in 
practice the degree to which radiological images can be restored to remove geometric 
blurring, using conventional Fourier deconvolution strategies, is severely limited by noise. 
The relatively poor signal-to-noise ratio of radiological images is a direct result of the 
requirement to operate within the ALARP principle, so that patient radiation doses are 
'kept as low as reasonably practicable consistent with the intended purpose' (JIRMER 
200(1 ) . Radiological images are therefore expected to be inherently noisy, with quantum 
noise being the dominant noise source. The noise power extends to the highest spatial 
frequencies recorded by the receptor, and it is the modulation at these higher spatial 
frequencies which is amplified using Fourier-based tech niques, result i ng in an image 



dominated by high frequency noise. The Wiener filter (" Wiener 



1949; 



HeLstronJ 11967) 



is a popular deconvolution tool which attempts to address the shortcomings of direct 
Fourier inversion techniques. Rs application to non-mammographic images has been 
demonstrated bv lDoughertv and KawaJ (|200l[ ). However, it is our experience that these 
techniques are of marginal value in radiology. 

In this work, an alternative image restoration approach based on a Bayesian 
maximum entropy method is applied to a radiological image deconvolution problem: 
reduction of geometric blurring in mammography. The aim of the work is to demonstrate 
an improvement in image spatial resolution with no associated penalty in terms of 
reduction in the signal-to-noise ratio perceived by the observer. The images used here 
were digitized ah initio] in principle digitized film could be used, but this would almost 
certainly be an unacceptable inconvenience if digital processing of this nature were 
to become part of clinical practice. Indeed it is the advent of Computed and Direct 
Radiography which we believe make this work particularly timely. 



2. Medical image restoration 

There is a large body of literature describing methods which aim to improve the quality 
of image information content through restoration techniques. Medical images are usually 
measurements of photon flux, and the data is usually noisy and often incomplete; 
statistical methods have yielded some robust methods of estimating the 'true' image 
distribution in these circumstances. The Expectation Maximization (EM) algorithm 
( Dempster et a\ 1977 ) has been widely applied to finding solutions which maximize a 
likelihood functi on for quantum noise limit e d data; in the me dical context this method 
has been used by Hudson and Larkin ( 1994 ) , De Pierro ( 1995 ) and Kinahan et al ( 1997 ) 
for tomographic ima ging. In the astronomical community EM is known a s Richardson- 
Lucy deconvolution ( Richardson 1972 : LucvHl974 : Shepp and Vardi 1982 ). Snvder et al 
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(1992) describe a method which maximizes the mean value of the log-likelihood for 
quantum noise limited data. It has been shown that this is equivalent to minimizing 
Csiszar's I-divergence (jCsiszarlll99ll ). a quantity equal to the negative of the entropy 
expression, —S, given in equation (J7J). The usefulness of Bayesian restoration stems 
from the fact that it allows the incorporation of sophisticated a priori knowledge into 
the formulation of the restoration method, while quite naturally enf orcing desirable 
properties such as positivity in the restoration. It has been argued by ISkiUinel (J1989) 
that in the absence of fur ther prior information entropy is the only consistent prior for 
positive, additive images. O' Sullivan et a\ (1998) gives a summary of these methods in 
terms of information theoretic image formation models. 



2.1. MEM theory 

The links betw e en st atistical mechanics and information theory were established by 
Javne? ( 1957q|Q l 1968h. I mage restoration using MEM was first described in a landmark 
paper by Friedenl ( 1972 ). The driving force behind practical implementations of the 
method came from radio astronomy and the nee d to improve radio maps of the sky 
f|Gull and Danielllll978i ICornwell and Evandll985h . 

MEM is a deconvolution technique derived from the forward map (1) for the imaging 
system, which relates postulated hidden data Xij to the observed data d^. 
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dij = * h + n i: j (1) 

Equation ([TJ is the forward map for the imaging system. The PSF, denoted by h, 
characterizes the geometric blurring of the imaging system, which acts on the hidden 
data. The hidden image is then further corrupted by additive noise to produce the 
observed image d^. 

The goal of image restoration is to obtain a solution x , which approximates to the 
hidden image x, as closely as the data and noise allow. We outlin e briefly the M EM 



approach to the problem here, but for a more complete review see iJavnesI (119821 ) and 
Skilling and Gull (jl985h . The notation in the following analysis has been simplified in 



relation to equation ([T|l. in that the data etc. are represented as vectors rather than 2D 
arrays. 

A trial restoration x is obtained and used as an initial guess for the hidden image 
x. The trial restoration is blurred, via equation (P), to generate mock data d. The x 
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goodness of fit statistic is used to measure the degree of misfit between the observed 
data and mock data: 

2 _ Yliidj — dj) , . 

X - - 2 W 

where a 2 is the variance in the noise, here taken to be constant, ft might be thought 
that a good approach would be to minimize this degree of misfit by choosing a suitable 
x, but this process is equivalent to the straightforward matrix inversion: 

x = A~ l d (3) 

where A is the matrix representation of the linear forward map equation ([!}. The 
problem here is that typically the matrix A is ill-conditioned, i.e. almost singular, so 
for a given d there are many vectors x which almost satisfy equation (JHJ), not necessarily 
close together. Thus if d contains even a small amount of noise, the resulting x can 
be far from the true image x; in other words the reconstruction is then dominated by 
noise rather than data, often referred to as over-restoration. In the presence of noise 
we would not in any case expect \ 2 to be minimized, but rather to be reduced to the 
appropriate value \ 2 = N, where N is the number of pixels in the image. Even when A 
is not ill-conditioned there are many x which satisfy this criterion. MEM is an example 
of regularized deconvolution and relies on a scheme of iterated forward maps rather than 
attempting to find a direct solution of the inverse problem. 

MEM treats the restoration process as a statistical inference problem based on 
Bayes' theorem and the aim is to obtain the most probable image x given the data: 

P(x | d) oc P(d | x) x P(x) (4) 

The likelihood P(d | x) is determined from our knowledge of the forward map equation 
the image noise is mainly quantum (photon) noise, which is modelled as a Gaussian 
process, so the likelihood term is quantified by the \ 2 distribution: 

P{d | x) oc exp(-x 2 / 2 ) (5) 

The prior P(x) is the probability that would be assigned by to a particular reconstruction 
x prior to the introduction of observational constraints. It is assumed that each 
luminance quantum has an equal a priori chance of falling into any pixel, in which 
case it is easy to show that 

P(x) oc exp(S) (6) 

where S is the configurational entrop y of the hidd en image. For positive additive 



distributions the entropy is defined as (jSkillind ll988): 



S{x) =^2[£i -m - Xi\n(xi/m)} (7) 

i 

Here m is a (constant) default level; S is maximized when £i = m, giving a flat featureless 
reconstruction. If each luminance quantum has a different chance of falling into any 



Mammographic image restoration using maximum entropy deconvolution 



5 



pixel, the default levels m; for pixel i can be chosen accordingly, and S is maximized 
when £i = rrti. A non-flat default level would be appropriate for example when the 
X-ray illumination is not uniform or when other knowledge of the restoration is known 
a priori. The strict definition of posterior probability is therefore: 



P(x | d) oc exp(S - x 2 / 2 ) 



(8) 



Values of £i which maximize this probability should be determined. This cannot 
be achieved analytically, and a numerical scheme must be employed. In practice a 
pragmatic modification of equation (JHJ) is used: 



P(x | d) oc exp(aS' — % 2 ) = exp(Q) 



(9) 



where Q = aS — x 2 , and a is a multiplier, as yet undetermined. x 2 is a measure of the 
misfit between the actual observations and those corresponding to the trial values xf, 
S is a measure of the structure within the trial image. If a is too small then too much 
weight will be given to the data, thus amplifying the effects of noise; if a is too large 
then too much weight will be given to the entropy and each Xj will be dragged towards 
its default value m^, thus losing real features. The procedure is to maximize Q with 
a initially fixed; a particular value of x 2 w iU correspond to this maximum. The aim 
is that this value should coincide with the expectation \ 2 = N\ the value of a which 
achieves this aim is found by a process of iteration. 

Initial feasibility tests were undertaken using the MATLAB Optimization Toolbox 
(The Mathworks Ltd, Cambridge, UK), particularly the constrained nonlinear 
minimization function fmincon in its LargeScale configuration, with a positivity 
constrain t on each £j. This function employs a preconditioned conjugate gradient 



method (jColeman and Lil Il994l Il99(l ). Our procedure has two loops. The inner 
one, in which a is fixed, iterates towards the minimum value of — Q for that a, 



and generates a corresponding value \ 2 



a 



The second loop iterates over a to 



minimize x 2 { a ) ~~ N , and is terminated when y 1 (a) falls within the narrow statistical 



range N ± (2iV) 1/2 jlackson and Ward! Il98 lr ) . Typically each loop requires 20 to 
30 iterations, and the final value of a is of order 20. After our successful initial 
implementation commercially available C++ software was purchased (MEMSYS5, 
Maximum Entropy Data Consultants, Bury St Edmunds, UK). This employs essent ially 
the same scheme, but is highly optimized, as described in S killing and Bryanl (jl984f ) and 
Gull and Skillind (J1999J). The MEMSYS5 software can handle larger images than our 
initial implementation, with much shorter processing times. 



3. Experimental Method 

Image quality comparisons were carried out using the Leeds TORMAM test object 
(figure C]) at various settings of geometric magnification. This phantom contains three 
groups of test features; fibres, simulated microcalcification clusters and low contrast 
plastic disks, plus an area designed to give an anthropomorphic impression of a breast 
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Figure 1. Layout of the quantitative side of the Leeds TORMAM test object 



parenchymal pattern with overlying microcalcification clusters (jCowen et a ll992h . This 
latter area of the phantom was not used in the present work. The phantom was imaged 
in various configurations on a laboratory-based Siemens Mammomat 3 mammography 
unit, with measured focal spot sizes of 0.7x0. 3mm broad focus and 0.2x0. 2mm fine 
focus. In all cases the phantom was imaged on top of a 3.5cm thick stack of Perspex 
to provide realistic scatter and attenuation as recommended in the instructions for the 
phantom. The radiographic factors used were 28kVp, 40mAs. These factors were chosen 
to be representative of the values used in routine mammographic quality assurance tests. 
The level of quantum noise in the test images was therefore realistic. The digital image 
receptor was part of a Philips ACR-3 Computed Radiography system (100 //m pixel 
size), comprising 1770 x 2370 pixels with overall dimensions 18 x 24 cm. 

In order to provide a sample of the point spread function (PSF) in the plane of 
the phantom, a piece of brass foil with a pinhole was included in the test images. In 
any practical application of this technique, the PSF information would be obtained from 
calibration images appropriately scaled for the position of the object. Including the PSF 
with the image was convenient in this demonstration as the scaling and interpolation 
step was avoided. The measured overall PSF is a convolution of one due to geometry 
and the the detector PSF, but is nevertheless the appropriate one to use. In fact the 
contribution from the detector PSF was small, evidenced by the fact that the measured 
PSF scaled with focal spot size and magnification in the appropriate manner (see figure 

E}. 

Three imaging geometries were used to illustrate varying degrees of focal spot 
geometrical blurring: 
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Magnification 1.8, broad focus (1.8BF). The plane of the phantom was 33.5cm 
from the tube focus and 27cm from the plane of the receptor. The image would be 
expected to be unacceptably degraded by geometric blurring. 

Magnification 1.8, fine focus (1.8FF). The phantom was in the same position as 
above, but imaged using the fine focal spot. This is the conventional magnification 
view provided on this mammography unit and used clinically. 

Magnification 3.0, fine focus (3. OFF). The plane of the phantom was 20.5cm from 
tube focus and 40cm from the plane of the receptor. This high magnification factor 
is not supported on standard mammography units because even using fine focus, 
the image would be unacceptably blurred. 

Images of the focal spots obtained in each configuration are shown in figure 121 



(a) 1 .8BF PSF (b) 1 .8FF PSF (c) 3.0FF PSF 




4 8 12 4 8 12 4 8 12 



Figure 2. Raw images of the focal spot obtained under the configurations (a) 1.8BF, 
(b) 1.8FF and (c) 3. OFF. 



The original and processed sets of images were viewed and scored by two 
independent observers both experienced in the use of mammographic image quality test 
phantoms. The test images were graded using the 3, 2, 1, scoring system recommended 
in the TORMAM phantom instructions and adop ted in surveys of mammogra phic image 
quality in the UK Breast Screening Programme ([Young and Ramsdale 2003 ). 



4. Image processing procedure 

The original images were in Siemens SPI file format which were 12-bit grayscale images 
holding 10-bit image data. A conversion to DICOM format was necessary to allow 
MATLAB to read the image data, but the conversion process retained all the luminance 
information held in the original files. 

MATLAB R12 was used to perform some preprocessing so that the subsequent 
MEM processing could be applied. These steps were: 
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Distance along TORMAM phantom 

Figure 3. 1.8BF. Pixel intensity profiles through particles in group A. The 
preprocessing step reduces and flattens the varying background of the image but leaves 
the structure and noise intact. MEM processing smoothes the noise and increases the 
signal-to-noise ratio of features deemed statistically to be caused by real objects. 



The original images have a varying background brightness caused by non-uniform 
intensity of the X-ray beam (which is by design in mammographic systems) and 
by variations in sensitivity over each digital plate (several such were used). To 
make the im ages more amenable to MEM processing it was desirable to reduce this 
background ( Donoho et a,\ 1992 ) . A Gaussian low pass filter was constructed and 
applied to the original images in the Fourier domain to obtain a background map. 
This map was subtracted from the original leaving a flattened image, which retained 
the phantom's high frequency features (see figure EJ). It should be noted that the 
Philips ACR-3 system is not a linear one, as each pixel value is proportional to 
the logarithm of the number of photons striking the detector; however, as the raw 
signal to mean background ratio is small (see figure EJ), typically less than 0.1, this 
procedure effectively linearizes the system, necessary for application of any of the 
methods discussed in section |21 

Removal of spurious, bright pixels not corresponding to any real feature in the 
TORMAM phantom. These pixels would lead to artefacts in the restored image. 
The removal was done by observing the image histogram and setting high valued 
outliers to some mean background level. 

The PSF images were cropped from the original TORMAM images. In each case 
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the noisy background of the PSF was filtered out (typically by discarding those 
pixels with less than 10% of the PSF peak value). This had the effect of slightly 
narrowing each PSF - thus leading to a conservative under-restoration. 

MEMSYS5 was used to treat the images following the MATLAB preprocessing. 
The interface to the MEMSYS5 kernel accepts image and PSF-image files as inputs 
(MATLAB image format was convenient) and allows the setting of certain parameters 
related to the theory described in section 2. The default image model m was defined to 
be a flat image with low pixel luminance values: our belief that, in the absence of data, 
very high photon counts were recorded (i.e. no absorption due to intervening material). 
MEMSYS5 uses a slightly more sophisticated forward map that that presented in 
equation ((TJ). 

The processed images contained between 2.2 million and 3.9 million pixels. 
MEMSYS5 typically converged to a solution within 15 a-iterations with a processing 
time of four to eight minutes, for an image of the full test object shown in figure 
The smaller images to be presented in figures |U to El were cropped from such an image 
after processing, rather than processed individually. The processed output files from 
MEMSYS5 were 8-bit PNG files. All image processing was carried out on a Pentium 4 
2.4GHz machine with 512MB of RAM. 



5. Results 



Before proceeding to a systematic comparison of the original and MEM processed 
images and related scores, we will say a few words about our initial aspirations for this 
technique. As mentioned in Section 1 these were to show improved spatial resolution 
without reduction in signal-to-noise ratio perceived by the observer, particularly a 
level of performance in which MEM processed 1.8BF images are at least as good as 
unprocessed 1.8FF ones. Our belief is that such a performance would be of clinical 
interest . Figures |U and El compare appropriate images of a filament group and simulated 
microcalcifications. These figures show that MEM reconstruction can give improvements 
in both resolution and perceived signal-to-noise ratio, and that the MEM processed 
1.8BF images are marginally better in this instance than the 1.8FF unprocessed ones in 
both respects. This promise encouraged us to undertake a systematic evaluation, and 
to consider magnifications greater than those normally used in clinical practice. Three 
particle groups (B,D and E) were chosen for this evaluation. The nominal particle 
size ranges for the se groups are 180-283 jum, 106-177 /im and 90-141 /im respectively 
(Co wen et a ll992h . In addition to the scoring system mentioned earlier we also calculate 
some signal to noise ratios (SNR) for these particular features. 

Figure El shows images of the three particle groups, cropped from the images taken 
under 1.8BF conditions. The selected particle group is clearly visible in figures El (a) 
and (d), discernible in figure El (e), just discernible in figure El (b) and not seen in figures 
El (c) and (f); however, we show the latter images to preserve the two by three format, 
which eases comparison with later improvements. As expected, unprocessed images 
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(a) 1 .8FF filament group C: Original (b) 1 .8BF filament groups C: MEM 




Figure 4. Comparison of the filament group C imaged under 1.8FF with the same 
feature imaged under 1.8BF and processed with MEM. 



(a) 1 .8FF particle group B: Original 



(b) 1 .8BF particle group B: MEM 





Figure 5. Comparison of particle group B imaged under 1.8FF with the same feature 
imaged under 1.8BF and processed with MEM. 



obtained in this configuration are too blurred to be clinically useful. However, the 
MEM processed images show significant improvements in resolution and signal-to-noise 
ratio, particularly regarding the high frequency noise, which were quantified using cuts 
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similar to that shown in figure EJ the quoted SNR is the difference between the largest 
signal and the mean background within each group, divided by the standard deviation 
of the signal in a region close to the group. Figures for Group B are SNR orig = 9.0 and 
SNRmem = 50.3. Group D: SNR orig = 5.9 and SNR M em = 13.7. Particle group E is not 
detected at this setting. 

Figure [7| shows cropped images of the same three particle groups, obtained with the 
conventional 1.8FF configuration. Remarks regarding group visibility are as for figure El 
The pinhole PSF imaged under these conditions is small and approximately Gaussian in 
shape, comprising just a few pixels. As expected, images obtained in this configuration 
are sharper than in the 1.8BF case; nevertheless MEM restoration still shows significant 
improvements. SNR measurements for Group B are SNR orig = 9.0 and SNR M em = 49.8. 
Group D: SNR orig = 5.4 and SNRmem = 18.5. Particle group E is not detected at this 
setting. 

Figure |H1 shows cropped images of the same three particle groups, obtained with 
an unconventional 3. OFF configuration. In this case all three groups are detectable 
in the original image of the phantom, but the MEM processed images show clear 
improvements in resolution, enabling fine details of individual microcalcifications to 
be discerned in image (d). SNR measurements for Group B are SNR orig = 8.8 and 
SNRmem = 49.8. Group D: SNR orig = 5.8 and SNR ME m = 22.4. Group E: SNR orig = 5.2 
and SNRmem = 10.4. 

To effect the same comparison as that illustrated in figure El between unprocessed 
1.8FF images and MEM processed 1.8BF ones, the bottom row in figure El should be 
compared with the top row in figure 

Figure El shows images of filament groups B, D and F taken with the 1.8FF 
configuration. These respectively have diameters 0.35, 0.25 and 0.20 mm, length 10 
mm. All three groups are detectable in the original and MEM processed images, but 
again the latter show improvements in resolution and visibility. To avoid information 
overload we do not show the 1.8BF and 3.00FF cases here, but experience with these 
matches that with the particle groups, and is quantified in the table of scores discussed 
below. 

The results obtained from scoring the original and MEM processed images under 
each of the three imaging geometries are shown in Table The processed images show 
an improved score in each case. This implies that the signal-to-noise ratio perceived by 
the observers for the various test features, including low contrast objects, was increased 
by the application of the MEM processing. 

6. Discussion 

The aim of the work was to demonstrate an improvement in spatial resolution for 
realistic radiological images from MEM de-blurring, with no associated penalty in 
terms of reduction in the signal-to-noise ratio perceived by the observer. The expected 
resolution improvements are shown in figures E1E1 Somewhat unexpected, however, 
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Figure 6. 1.8BF. Close-ups of particle groups B, D and E are shown in (a), (b) and (c) 
respectively. MEM restorations are shown beneath in (d), (e) and (f); improvements 
in visibility and resolution are apparent with the first two groups. Group E is not 
detectable in either the original or MEM processed images. 
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Figure 7. 1.8FF. Close-ups of particle groups B, D and E are shown in (a), (b) and (c) 
respectively. MEM restorations are shown beneath in (d), (c) and (f); improvements 
in visibility and resolution are apparent with the first two groups. Group E is not 
detectable in either the original or MEM processed images. 
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Figure 8. 3. OFF. Close-ups of particle groups B, D and E are shown in (a), (b) and (c) 
respectively. MEM restorations are shown beneath in (d), (e) and (f); improvements in 
visibility and resolution are apparent in all three cases, with fine details of individual 
microcalcifications being discernable after MEM restoration (d). 




Figure 9. Close-ups of filament groups B, D and F are shown in (a), (b) and (c) 
respectively. MEM restorations are shown beneath in (d), (e) and (f); improvements 
in visibility and resolution are apparent in all three cases. 
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Table 1. Image scores of the original and MEM processed images, obtained by 
averaging the individual scores of the two independent observers. 



Mode 


Image 




TORMAM 


scores 




Filaments 


Particles 


Disks 


Total 


1.8BF 


Original 


20.5 


6.0 


24.5 


51.0 




MEM 


38.5 


8.5 


30.0 


77.0 


1.8FF 


Original 


27.0 


6.0 


21.5 


54.5 




MEM 


41.0 


10.5 


29.5 


81.0 


3.0FF 


Original 


42.0 


9.5 


27.5 


79.0 




MEM 


57.5 


11.0 


35.5 


104.0 



were the improvements in image score shown in Table 1, as most of the features in 
the TORMAM phantom, i.e. the fibre groups and disks, are essentially low-contrast 
features whose detection would be expected to be limited by the relative noise level in 
the image. The improvement in scores therefore implies an improvement in signal-to- 
noise ratio for this phantom. For the filaments and particles, which are comparable 
to the PSF in extent, improvements in visibility are effected by enhanced intensity and 
sharpness due to focussing, and by noise reduction. For the disks, which are significantly 
larger than the PSF, sharpness (resulting in easier edge detection) and noise reduction 
are the important factors. 

It may be that the improvements in features imaged against a uniform background, 
as in these demonstrations, are better than those which might be achieved when imaging 
diagnostic features against an anatomical background. Suitable experiments are under 
way. 

There are a number of possible applications of MEM deconvolution in radiology. 
This paper has demonstrated the use of high geometric magnification in conjunction 
with a relatively large focal spot. Other options might be to use larger focal spot tubes 
to increase heat capacity and therefore allow extremely short exposures, or high outputs 
which could be used with high filtration. Other deconvolution problems in radiology 
include light scatter in image intensifier optics and reduction of the effects of scattered 
radiation. 

7. Conclusion 

An image processing approach based on a maximum entropy method has been applied to 
the problem of restoring focal spot geometric blurring in magnification mammography. 
The results show an improvement in image spatial resolution and an improvement in 
terms of the image signal-to-noise ratio perceived by the observer, as evaluated using a 
standard phantom and at a realistic quantum noise level. This successful demonstration 
of image de-blurring in noisy radiological images offers the possibility of weakening the 
link between focal spot size and geometric blurring in radiology, and thus opening up 
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new approaches to system optimization. 
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